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Abstract. Sub-Saharan Africa harbours an outstanding diversity of tortoises of which the leopard tortoise Stigmochelys 
pardalis is the most widespread. Across its’ range the species is impacted by habitat transformation, over-collection for hu- 
man consumption and the pet trade, road mortality, and electrocution by electric fences. Most leopard tortoises in south- 
ern Africa are nowadays restricted to reserves and private farms. So far confiscated tortoises are frequently released into a 
nearby reserve without knowledge on their area of origin. This is problematic, as it has been demonstrated that the leopard 
tortoise harbours five distinct mitochondrial lineages, of which three occur in the southern portion of the species’ distri- 
butional range (South Africa, Namibia, and Botswana). Using 14 microsatellite loci corresponding to 270 samples collect- 
ed throughout southern Africa, we found a clear substructuring in the north constituting four clusters (western, central, 
north-eastern, and eastern). Genetic diversity was particularly high in the north-east and decreased towards the south. In 
addition, we found a significant size difference between the studied populations. Our basic morphological analysis showed 
that tortoises from the southern cluster tend to grow bigger than tortoises from the north. We established a comprehensive 
genetic database for South Africa and Namibia that can serve as a conservation management tool for the assignment and 
potential release of translocated or seized leopard tortoises based on genetic affiliation. 
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Introduction 


Tortoises and turtles are among the most imperilled bio- 
ta on the planet, with approximately 60.4% of 356 species 
threatened or recently extinct (Lovicu et al. 2018, Turtle 
Conservation Coalition 2018). Of 121 tortoises that ex- 
isted since the Pleistocene, 69 (57%) already disappeared 
(Lovicu et al. 2018). The impacts that are responsible for 
this alarming situation include habitat destruction, cli- 
mate change, unsustainable overexploitation for human 
consumption as well as for the international pet trade, and 
the introduction of pathogens (WIMBERGER et al. 2011, 
Lovicu et al. 2018). These also apply to the leopard tortoise 
Stigmochelys pardalis (BELL, 1828) (EVANS 1988). Between 
1987 and 1991, S. pardalis accounted for 76% of tortoises 
exported from Africa for the international pet trade and 
declines in some areas have been attributed to unsustain- 
able collection and trade (BRANCH 2012). Over the last five 
years 171,444 living S. pardalis, have been exported from 
Africa for commercial purposes (CITES trade database 
2020). As a consequence, the once widespread species is 
now considered threatened in some regions, while it has 


already disappeared from few areas in eastern and west- 
ern South Africa (HOFMEYR & BAARD 2014). Natural bush 
fires and deliberate fires set to promote regrowth and in- 
crease forage quality for livestock account for a major pro- 
portion of injuries and mortalities (KABUGUMILA 2001), 
while severe droughts and road mortality also contribute 
to population decline (BOYCOTT & BOURQUIN 2000). To- 
day the vast majority of tortoises, including S. pardalis, are 
found on privately owned fenced farm land and reserves. 
Unfortunately, their large body size and highly domed car- 
apace make leopard tortoises susceptible to being killed by 
electric fences (BURGER & BRANCH 1994, ARNOT & MOLTENO 
2017, MACRAY 2017). 

Recent molecular genetic analyses revealed the species to 
constitute five distinct mitochondrial lineages (SPITZWEG et 
al. 2019). One is distributed in the south, a second in the 
north-west, and a third in the north-east of southern Africa, 
while two additional lineages are distributed further north 
(SPITZWEG et al. 2019). Using microsatellites SprrzweEe et al. 
(2019) revealed the southern mitochondrial clade to match 
with a well-defined southern nuclear cluster, whilst the 
north-western and north-eastern clades of southern Africa 
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corresponded to another nuclear cluster with three subclus- 
ters. Besides the genetic diversity, the leopard tortoise ex- 
presses considerable geographic variation in body size, shell 
shape, and colouration. LovERIDGE and WILLIAMS (1957) 
attributed these differences to subspecies variation. How- 
ever, the distribution of the mitochondrial lineages or the 
respective nuclear clusters does not match with the recog- 
nized subspecies (FRITZ et al. 2010, SPITZWEG et al. 2019). 
These findings demonstrate the importance of the applica- 
tion of conservation genetic techniques to preserve the ge- 
netic diversity of the species. While previous studies suf- 
fered from incomplete sampling in some areas (SPITZWEG et 
al. 2019) our present study adds crucial data from northern 
South Africa completing these sampling gaps. 

We examined the population structure and reveal ge- 
netic differentiation of S. pardalis using fine scale analyses 
of nuclear microsatellite markers. We deliver an essential 
database for future conservation relocation practices that 
will improve the conservation management and preserva- 
tion of the genetic diversity of S. pardalis. 


Material and methods 
Sampling, selection of microsatellite loci, 
and general data evaluation strategy 


In total, 270 leopard tortoises from all across South Africa, 
adjacent Namibia, and southern Mozambique were studied 
(Supplementary document S1). These include 204 speci- 
mens previously sampled (SprrzwEc et al. 2019). For genet- 
ic analyses, a small amount of blood or tissue was collected 
from live tortoises in accordance with methods approved 
by the Ethics Committee of the University of the Western 
Cape (Ethics Reference Number (AR 19/4/1). All speci- 
mens were subsequently released at their respective col- 
lection site. Blood samples were preserved on FTA classic 
cards (Whatman, GE Healthcare, Munich, Germany) and 
stored at room temperature. For most tortoises sampled in 
South Africa, the carapace length (CL) was measured to 
the nearest 1 mm using a measuring tape. Fieldwork and 
sampling in South Africa was permitted by the Limpopo 
Provincial Government (ZA/LP/91608), the Department of 
Environmental Affairs, Biodiversity Northern Cape Prov- 
ince (245/2015), Ezemvelo KZN Wildlife (OP 139/2017), and 
Biodiversity North West Province (NW 6124/10/2018). 

The same 14 nuclear microsatellite loci were targeted as 
in our previous study on S. pardalis (SPITZWEG et al. 2019). 
For DNA extraction, PCR, and genotyping information see 
SPITZWEG et al. 2019. 


Data analyses 


For conservation genetic purposes microsatellite data 
of 270 leopard tortoises were analysed using the unsu- 
pervised Bayesian clustering approach of STRUCTURE 
2.3.4 (PRITCHARD et al. 2000, Husisz et al. 2009), apply- 
ing the admixture model and correlated allele frequencies. 
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STRUCTURE searches for partitions which are, as far as 
possible, in Hardy-Weinberg equilibrium and linkage equi- 
librium. Unsupervised analyses were selected because this 
approach clusters samples strictly according to their genetic 
information, without making presumptions about popula- 
tion structuring (e.g. geographic distances, sampling sites). 
All calculations were performed for K = 1-10, and the most 
likely number of clusters (K) was determined using the AK 
method (EvANNO et al. 2005) as implemented in the soft- 
ware STRUCTURE HARVESTER (EARL & VON HOLDT 
2012) and mean Ln probabilities. Calculations were repeat- 
ed 10 times for each K using a MCMC chain of 750,000 
generations and a burn-in of 250,000 generations. Popu- 
lation structuring and individual admixture were visual- 
ized using DISTRUCT 1.1 (ROSENBERG 2004). Individuals 
with a membership proportion below 80% were considered 
to have admixed ancestries (BARILANI et al. 2007, RANDI 
2008). STRUCTURE is prone to bias from uneven sample 
sizes (PUECHMAILLE 2016) and typically detects only the 
uppermost hierarchical level of population differentiation 
(EVANNO et al. 2005). Therefore, the analysis was repeated 
for subsamples corresponding to previously identified clus- 
ters but excluding admixed individuals (Supplementary 
document S1). In addition, STRUCTURE was repeated for 
each of the respective subclusters of the northern cluster. 

Further, Principal Component Analyses (PCAs) were 
performed to examine population structuring using mi- 
crosatellite data and the package ADEGENET (JOMBART 
2008) for Cran R 3.2.3 (R Development Core Team 2015). 
PCAs are less sensitive to sample size and are independ- 
ent from population genetic assumptions (PUECHMAILLE 
2016). Two distinct PCAs were executed. In the first anal- 
ysis, all 270 individuals were included, and symbols were 
coloured after STRUCTURE clusters (north, south, ad- 
mixed). The second PCA included only the tortoises from 
the northern STRUCTURE cluster (151 individuals), with 
symbol colours corresponding to STRUCTURE subclus- 
ters (N1, N2, N3, N4, Admixed). 


Diversity within and divergence 
between population clusters 


Diversity and divergence were estimated for population clus- 
ters revealed by STRUCTURE. A frequency table for mic- 
rosatellite alleles was produced using COVERT 1.31 (GLAv- 
BITZ 2004). ARLEQUIN 3.5.2.1 (EXCOFFIER & LISCHER 
2010) was applied to infer locus specific observed (H,) and 
expected (H,) heterozygosity and to perform AMOVAs 
(10,000 permutations) for microsatellite clusters revealed 
by STRUCTURE. FSTAT 2.9.3.2 (GOUDET 1995) was used 
to compute values for locus specific allelic richness. 


Morphology 


Data on carapace length (CL) of 137 adult tortoises was an- 
alysed using the packages DPLYR (WIckHAM et al. 2019), 
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TIDYR (WICKHAM & HENRY 2019), and MULTCOMP 
(HoTHoRN et al. 2008) and visualized using the packag- 
es GGPLOT2 (WICKHAM 2016), GGTHEMES (ARNOLD 
2019), and KNITR (X1z 2019) for Cran R 3.2.3. A two-way 
ANOVA was applied using main genetic clusters identified 
by STRUCTURE (south/north) and sex as fixed effects to 
test for significant difference in carapace length (response 
variable). Subsequently, differences were identified using a 
Tukey HSD test (DuBITZKY et al. 2013). 


Results 
Population structuring 


Using all samples (N = 270) STRUCTURE revealed two 
clusters as the best solution (Supplementary documents S2, 
S3), one corresponding to southern sampling sites and the 
other to more northern sites within South Africa, adjacent 
Namibia, and southern Mozambique. Admixed ancestry 
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was only observed in a few sites in the Provinces Western 
Cape, Northern Cape, and North West (grey coloured sec- 
tion of pie charts; Supplementary document S2). In addi- 
tion, tortoises from both clusters were found to live in close 
proximity to each other in several localities (Supplemen- 
tary document S2). 

Additional STRUCTURE runs for each of the two clus- 
ters excluding admixed individuals revealed K = 1 as the best 
solution for the southern cluster but suggested further struc- 
turing for the northern subset (K = 4; Fig. 1, Supplementary 
document S1). Thus, the northern subset (Fig. 1b: N1) includ- 
ed a cluster corresponding to samples from Namibia and 
the South African Northern Cape Province. Two additional 
clusters (Fig. 1b: N2, N4) contained samples from central- 
southern Africa with a clear north-south separation. The 
fourth cluster (Fig. 1b: N3) corresponded to a smaller area in 
the north-east of South Africa, and southern Mozambique. 
A high number of tortoises from central and north-eastern 
South Africa showed admixed ancestries (Fig. 1b). 


À 0 200 400 
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Figure 1. a) Genotypic clustering of the two subsets (northern and southern, excluding admixed individuals) as inferred by STRUC- 
TURE. The runs with the best probability values are shown. Distinct clusters are colour coded corresponding to STRUCTURE inferred 
clusters (see b). Within each cluster an individual is represented by a vertical bar that reflects its ancestry. Mixed ancestry is indicated 
by differently coloured segments corresponding to inferred genetic percentages of the respective cluster. b) The distribution range 
of S. pardalis is displayed as shaded area according to BAKER et al. 2015 and Sprrzwec et al. 2019. Individual sampling sites were 
lumped for clarity. Symbol sizes refer to the respective sample size. Slices represent percentages of individuals with mixed ancestries 
or conflicting cluster assignments. For original sites see Supplementary document S1. Inset shows the whole putative distributional 


range (black) of Stigmochelys pardalis according to TT WG (2017). 
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STRUCTURE runs for each of the four clusters of the 
northern subset revealed K = 1 as the best solution in all 
four runs (Supplementary document S4). Therefore, no ad- 
ditional substructure could be found in the north. 

Principal Component Analyses (PCAs) supported the 
results of the STRUCTURE runs and the near absence of 
admixture between tortoises of the northern and southern 
cluster, but indicated a high level of admixture between the 
four clusters of the northern subset (Fig. 2). 


Diversity within and divergence 
between STRUCTURE Clusters 


In agreement with its nuclear genomic substructuring 
(Supplementary document S1, 2), genetic diversity was 
higher among tortoises from the northern cluster than it 
was in the south (Table 1). However, despite similar sample 
sizes, the northern cluster had a higher haplotype diver- 
sity. For microsatellites the average number of alleles per 
locus (n, = 18.6) in the north was twice as high as for the 
southern cluster (n, = 9.3). The northern cluster had al- 
most twenty times more private alleles (n, = 138) compared 
to the southern cluster (n, = 7). The observed pattern was 
also reflected by a higher heterozygosity in the north (H,= 
0.709) than in the south (H, = 0.549). 

When comparing the four clusters of the northern sub- 
set (N1-N4, Fig. 1) the highest genetic diversity was ob- 
served within subcluster N2. With an average of 10.8, sub- 
cluster N2 had the highest number of alleles per locus, fol- 
lowed by subcluster N4 (9.6), Ni (8.4), and N3 (6.8). In 
concordance, subcluster N2 also had the highest number 
of private alleles (n, = 28) while subcluster N3 had the low- 
est (n, = 8). 

The two main clusters (northern and southern, Sup- 
plementary document S2) differed by an F,, value of 0.145 
(p < 0.05) while the northern subclusters N1-N4 (Fig. 1) dif- 
fered by substantially lower F„„ values ranging from 0.039 


Table 1. Genetic diversity of STRUCTURE clusters excluding 
admixed individuals. n - number of individuals; n, - number 
of alleles; n a -average number of alleles per locus; n, — number 
of private alleles; AR - allelic richness; H, - average observed 
heterozygosity; H, - average expected heterozygosity. 


Cluster n n, ny n, AR Ho H, 

North 151 130 186 138 17.611 0.709 0.803 
South 110 261 9.3 7 9.237 0.549 0.622 
N1 19- , 117 8.4 17 7.245 0.748 0.784 
N2 15 151 10.8 28 9.954 0.773 0.841 
N3 12 95 6.8 8 6.786 0.685 0.721 
N4 28 134 9.6 22 9.297 0.675 0.743 


to 0.107 (Table 2). In accordance, an AMOVA for the north- 
ern and southern clusters revealed 14.49% of the molecular 
variance to occur among and 85.51% within the two clusters. 
Regarding the northern subclusters Nı-N4 (Fig. 1) only 
6.24% of the molecular variance was found to occur among, 
while 93.76% of variance occurred within subclusters. 


Morphology 


A two-way ANOVA was conducted that examined the ef- 
fect of genetic cluster and sex on body size of 137 S. parda- 
lis. The northern cluster was composed of 42 males (mean 
CL = 34.3 + 7.9 cm SD) and 28 females (mean CL = 40.8 + 
9.3 cm SD) whereas the southern cluster contained 35 
males (mean CL = 45.4 + 7.3 cm SD) and 32 females (mean 
CL 42.7 + 10.9 cm SD). We found a statistically significant 
difference in carapace length between clusters (n unen = 67 
Donnen = 703 two-way ANOVA: F (1) = 22.79, p = > 0.001) 
but not in sexes (naes = 7% Demas = 60; two-way ANO- 
VA: F (1) = 1.50, p = 0.236). A subsequent Tukey HSD post 
hoc test based on the results of the ANOVA interaction be- 
tween location and sex (F (1) = 9.04, p = > 0.01) revealed 


(A) (B) 
All North 
. 3 4 i 
fm 
Š Hl ne 
HS North Gl ns 


i 3 «8 GS South 
‘ z “es axis | J Admixed : i 
‘las tel 7" ei A a 


MN 


bi M 


Figure 2. Principal Component Analyses (PCAs) using microsatellite data of A) 270 individuals (all samples) and B) 151 individuals 
of the northern cluster coloured according to STRUCTURE results. For A) axes 1-3 explain 10.25%, 5.22%, and 4.47% of variance; 
for B) axes 1-3 explain 6.45%, 5.27% and 5.01%, respectively. Ovals indicate 95% confidence intervals. 
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Table 2. Fixation indices (F, values) for microsatellite data of the 
northern subclusters N1-N4 (Fig. 2). All values are statistically 
significant (p < 0.05). 


N1 N2 N3 


N2 0.039 - 
N3 0.093 0.051 - 
N4 0.059 0.031 0.107 


that male tortoises from the southern cluster were signifi- 
cantly larger (Tukey HSD test, p < 0.001) than male tor- 
toises from the northern cluster (Fig. 3). In contrast no 
significant difference was observed between northern and 
southern female tortoises (Tukey HSD test p = 0.84). With- 
in the northern cluster, females were found to be signifi- 
cantly larger (Tukey HSD test p = 0.02) than males while 
no intersexual difference could be demonstrated between 
females and males from the southern cluster (Tukey HSD 


test p = 0.59). 


Discussion 


In recent years wildlife forensics became a key means of 
enforcing legislation surrounding the illegal trade in pro- 
tected and threatened species (Mucci et al. 2014, OGDEN & 
LINACRE 2015), which also holds true in South Africa. How- 
ever, the identification of the geographic origin of seized 
animals is based on our ability to assign a sample to a par- 
ticular population. This requires the source population to 
be sufficiently genetically distinct from other populations 
and is usually reliant on the existence of population data 
from multiple geographic areas (OGDEN & LINACRE 2015). 
Our study provides a comprehensive genetic database for 
wildlife forensics of translocated or seized leopard tortois- 
es in southern Africa developed in collaboration with the 
South African National Biodiversity Institute (SANBI) and 
the Gauteng Department of Agriculture and Rural Devel- 
opment (GDARD). Compared to SPITZWEG et al. (2019), 
our sampling covers the whole distribution range of S. par- 
dalis in South Africa, closing previous sampling gaps in 
central and northern South Africa. Here, we confirm the 
pronounced north-south genetic differentiation in south- 
ern African leopard tortoises reported by Sprrzwec et al. 
(2019) (Supplementary document S2) and demonstrate the 
distribution of the southern cluster to stretch further to the 
north than anticipated. According to BAKER et al. (2015) 
there is a major distribution gap for leopard tortoises in the 
western part of the Northern Cape Province and the north- 
western part of the Western Cape Province. However, we 
confirmed S. pardalis from the southern parts of the North- 
ern Cape Province, well within the proposed distributional 
gap, and reveal these tortoises to belong to the southern 
cluster (Supplementary document S2). Even though many 
herpetologists believe that S. pardalis has been introduced 
in western South Africa (HOFMEYR & BAARD 2014), our re- 
sults and field observations (confirming viable population) 


support a natural occurrence in this region. In addition, 
we investigated the geographic variation in body size due 
to the observations of LovERIDGE and WILLIAMS (1957). 
We demonstrate a pronounced geographic variation in 
carapace length of male tortoises between the two major 
microsatellite clusters from the north and south (Fig. 3). 
Males of the southern cluster reach carapace lengths of up 
to 65 cm (mean CL = 45.4 + 7.3 cm SD) while males from 
the northern cluster stay significantly smaller (max CL = 
50 cm, mean = 34.3 + 7.9 cm SD; Fig. 3). Although S. par- 
dalis females generally grow larger than males (BRANCH 
2012) we could not confirm a significant difference in cara- 
pace length between southern and northern females. Most 
probably a combination of environmental factors such as 
aridity, seasonality, and food availability are responsible for 
shaping the morphological variation in S. pardalis. Simi- 
larly, LOEHR et al. (2007) reported aridification to be cor- 
related with lower average and maximum body size in fe- 
males of the world’s smallest tortoise, Chersobius signatus. 
It is important to determine the geographic origin of 
seized leopard tortoises, as they can be locally adapted to 
a specific climate and might not survive non-analogous en- 
vironmental conditions as demonstrated by WIMBERGER 
et al. (2009). Our basic morphological analyses show that 
body size is not necessarily a proxy for geographic origin, 
even though considerable geographic variation was expect- 
ed (LOVERIDGE & WILLIAMS 1957). Furthermore, releasing 
animals from different genetic clusters in inappropriate re- 
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Figure 3. Differences in carapace length of females and males 
from the major microsatellite clusters (south and north). Box- 
plots show median, minimum, and maximum. Violin plots show 
the probability density of data. 
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gions can result in extinction or genetic pollution of local 
gene pools (BUTLER 1994, FITZPATRICK & SHAFFER 2007). 
A similar case was suspected in France where certain popu- 
lations of the European pond turtle Emys orbicularis appear 
to be polluted with genes from Italian populations (VamM- 
BERGER et al. 2015), but the impact has not been investi- 
gated. While our results are in agreement with SPITZWEG et 
al. (2019) in finding no substructuring in the southern clus- 
ter (Fig. 1), our increased sampling revealed the northern 
cluster to harbour an additional fourth subcluster (N4) in- 
stead of three as previously assumed (SPITZWEG et al. 2019). 
The additional subcluster occurs in the eastern part of the 
Northern Cape, throughout the Free State, and in the south- 
ern part of the North West Province of South Africa (Fig. 1). 
A high number of tortoises sampled in central and north- 
eastern South Africa show admixed ancestries (Fig. 1; grey 
coloured sections of pie charts). The reason for this may be 
the close proximity to the O. R. Tambo International Air- 
port in Johannesburg, where tortoises are frequently confis- 
cated and released nearby. The phylogeographic pattern of 
S. pardalis parallels with that of the puff adder Bitis arietans 
(BARLOW et al. 2013) and freshwater terrapins of the genus 
Pelomedusa (VAMBERGER et al. 2018) supporting a natu- 
ral distribution pattern. The distribution of the southern 
cluster of S. pardalis corresponds well with subclade Ic of 
Pelomedusa galeata, while the northern subclusters N3 and 
N4 correspond to Pelomedusa subclades Ib and Ia (VAM- 
BERGER et al. 2018). The northern subclusters Ni and N2 
mirror the distribution of Pelomedusa subrufa (VAMBER- 
GER et al. 2018). Stigmochelys pardalis is one component of 
the South African biota whose phylogeographic structure 
was shaped by physiography and climate, as has also been 
demonstrated for Pelomedusa (VAMBERGER et al. 2018). To 
preserve the genetic diversity and geographic structure of 
S. pardalis, we propose five management units for southern 
Africa. These correspond to the two major microsatellite 
clusters and the northern subclusters N1-N4 (Fig. 1). Using 
our genetic database, the confiscated tortoises should be ge- 
netically tested and assigned to one of the five microsatellite 
clusters before being considered for reintroduction into the 
wild. If assignment to a genetic cluster is successful, release 
in a safe location within the distributional boundaries of 
the correct genetic cluster can be considered. However, a 
suitable release protocol should include a thorough pre-re- 
lease health screening to prevent the introduction of patho- 
gens into wild populations (WIMBERGER et al. 2011). This 
study represents the foundation for further conservation 
management decisions and future conservation actions for 
leopard tortoises in southern Africa. 
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